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Abstract 

In this paper I shall resume the results of the MACRO experiment on atmospheric muon neutrino 
oscillations. 


1 Atmospheric neutrinos and MACRO 

A high energy primary cosmic ray interacts in the upper atmosphere producing a large number of pions 
and kaons, which decay yielding muons and muon neutrinos; the muons decay generating muon and 
electron neutrinos. These neutrinos are produced in a spherical shell at 10 — 20 km above ground and 
travel towards the Earth. The can be detected by an underground experiment as muons, after their 
charged current interactions with the material inside or outside the detector {vp.N —> -|- ...). 

MACRO was a large area multipurpose underground experiment optimized to search for rare events 
in the cosmic radiation. The detector was located in the Gran Sasso laboratory, which is well shielded 
from downgoing cosmic ray muons by a mean rock thickness of 3700 m.w.e.. MACRO was composed of 
three sub-detectors: liquid scintillation counters for energy and T.O.F. measurements, limited streamer 
tubes for particle tracking and nuclear track detectors for magnetic monopoles and nuclearites searches. 
It had a modular structure: it was divided into six sections referred to as supermodules, each of which 
had a size of 12.6 x 12 x 9.3 m^. A cross section of the detector is shown in Fig. 1 on the left, with the 
different event topologies of interactions in or around the detector: Internal Upgoing (lU), Upgoing 
Stopping (UGS) and Internal Downgoing (ID), and Up throughgoing muons. Data were taken from March 
1989 till April 1994 with the detector under construction and from the middle of 1994 till the end of 2000 
with the full detector. 


2 Neutrino oscillations analysis 

Neutrino oscillations. In the simple case of only two neutrinos with flavour eigenstates Vr) the 
survival probability of a pure beam is: 

= P{v^ ^ i/^) = 1 - sin^ 2023 sin^ ^ 1-27 Ato — 

where Pii^^ —> Vt) is the oscillation probability, 023 is the mixing angle, (z/ 2 , 1 ^ 3 ) are the mass eigenstates, 
Am^ = m§ — m 2 (eV^), L (km) is the distance travelled by the neutrino from the production point to 
the detection point, and Ei, (GeV) the neutrino energy. 
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Figure 1: Left: Cross section of the MACRO detector and sketch of different event topologies induced by ^ /r 
charged current interactions. The black points and the black rectangles represent streamer tubes and scintillator 
hits, respectively. Right: Zenith angle distribution of up throughgoing muons (black points). The dashed line is 
the expectation for no oscillations (with a 17 % scale uncertainty band). The solid line is the fit for an oscillated 
muon flux which yielded maximum mixing and Am^ = 2.5 • 10~® eV^ [2]. 


Up throughgoing muons. The up throughgoing muons come from interactions below the 
detector. These neutrinos have an average energy of ~ 50 GeV M- The data (6.16 years of livetime) 
deviate in absolute value and in zenith angle distribution from Monte Carlo (MC) predictions without 
oscillations, see Fig. 1 on the right. Fig. 2 upper left shows the allowed region for Vr oscillations 

in the sin^20 — Am^ plane, computed according to for the up throughgoing muon events; our allowed 
region is compared with those obtained by the SuperKamiokande (SK) ^ and Soudan 2 Q experiments. 

Matter effects due to the difference between the weak interaction effective potential for muon neutrinos 
with respect to sterile neutrinos would produce a different total number and a different zenith distribution 
of up throughgoing muons |^. In Fig. 2 on the upper right the measured ratio between the nearly vertical 
events with — 1 < cosQ < —0.7 and the horizontal events with —0.4 < cosQ < 0 is shown as a black point. 
In this measurement most of the theoretical uncertainties on neutrino flux and cross section cancel. The 
remaining theoretical error combined with the experimental error is estimated to be 7%. We measured 305 
events with —1 < cosQ < —0.7 and 206 with —0.4 < cosQ < 0; the ratio is R = 1.48 ± 0.13stat ± O.lOsys- 
For Am^ = 2.5 • 10“^eU^ and maximal mixing, the expected value of the ratio for —> i/r is Rr = 1.72 
while for —> iZg > is Rsteriie = 2.16. —> Vs oscillations are disfavoured at 99% c.l. compared to the 

I'r channel 0, (|]. 

The oscillation probability is a function of the ratio Lj see Eq. 1. In MACRO the neutrino energy 
was estimated for each event by measuring the muon energy via Multiple Coulomb Scattering (MCS) 
of the induced muon in the horizontal absorbers in the lower MACRO (Fig. 1, left). The r.m.s. of the 
lateral displacement for a muon crossing the apparatus on the vertical is <tmcs — 10 cm / (GeV). 
A first analysis was made studying the deflection of up throughgoing muons with the streamer tubes 
in “digital mode”. This method could reach a spatial resolution of 1 cm which implies a maximum 
measurable energy of 10 GeV ||]. To improve the spatial resolution of the detector, a second analysis 
was performed with the streamer tubes used in “drift mode” and the resolution achieved was ~ 3 mm. 
The method allowed the separation of the up throughgoing muons in 4 subsamples corresponding to 
average neutrino energies of 12, 20, 50 and 102 GeV, respectively. The comparison of the 4 zenith angle 
distributions with the predictions of the no oscillations MC shows a disagreement at low energies, while the 
agreement is restored at the higher neutrino energies. The distribution of the ratios R = {Data/MCno osc) 
obtained by this analysis is plotted in Fig. 2 in the lower part as a function of logio{L/E,j) [Q. 


2 

















10 -'' 

„ 10-2 
CM 
> 

3^ 

10-3 

< 


10 '"^ 


10-5 


r 



MACRO,' 

yCSK 


SOUDAN 2^.^ 

' "s. 

_L 

_^_1_^_ 


0 0.2 0.4 0.6 0.8 1.0 

sin2 (29) 




Figure 2: Upper left: 90 % c.l. allowed regions for ^ Vt oscillations from np thronghgoing muons and 
comparison with the Sondan 2 and SK allowed regions. Upper right: ratio of events with — 1 < cos9 < —0.7 
to events with —0.4 < cos6 < 0 vs Am^s for maximal mixing. The black point is the measnred valne, the solid 
line is the prediction for ^ v-r oscillations, the dash-dotted line is the prediction for —> VstariU oscillations. 
Bottom: Data/MC vs L/Ei, for np throngoing muons (black circles) and for semicontained up-^ (black square). 
The muon energy was estimated by MCS and by MC methods. The shaded region represents the uncertainty 
in the MC prediction assuming sin^ 29 = 1 and Am^ = 0.0025 eV^. The horizontal dashed line at Data/MC=l 
is the expectation for no oscillations. 


Low energy muons. An analysis was done also for the lower energy samples. These low energy 
muons come from neutrinos which have an average energy of ~ 4 GeV Q. The Internal Upgoing (lU) 
muons come from interactions in the lower apparatus; for these events two scintillation counters are 
intercepted; the T.O.F. is applied to identify upward going muons. 

The upgoing stopping muons (UGS) are due to external interactions yielding upgoing muons 
stopping in the detector. The internal downgoing muons (ID) are due to z/^-induced downgoing muon 
tracks with vertex in the lower MAGRO, see Fig. 1 on the left. A summary table of the number of 
events is presented with MC predictions in Table |^. The data show a uniform deficit for the whole 
angular distribution with respect to predictions, ~ 50% for lU, 75% for ID + UGS; there is a good 
agreement with the predictions based on neutrino oscillations with the parameters obtained from the up 
throughgoing muons. The average value of the double ratio R = {Data/MC)iu/{Data/MC)iD+UGS 
over the measured zenith angle range is i? ~ 0.77 ± 0.07; the error includes statistical and theoretical 
uncertainties; i? = 1 is expected in case of no oscillations. 

In conclusion, MACRO atmospheric neutrino flux analysis gives strong evidence to the —> Vr 

oscillation hypothesis, with Am^ = 2.5 x 10“^ eV^ and maximal mixing. 
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Events 

MC-No oscillations 

R = (Data/MCnooBc) 

Up throughgoing 

809 

1122 ± 191 

(0.721 ± OmGstat ± 0.043,y. ± 0.123t^) 

Internal Up 

154 

285 ± 28sy. ± 71th 

(0.54 ± 0.04.tat ± O.OS.ys ± 0.13th) 

Up Stop + In Down 

262 

375 ± 37sys ± 94t^ 

(0.70 ± 0.04,tat ± 0.07,j;s ± 0.17th) 


Table 1: Summary of MACRO r'/i —> /i events in — 1 < cos9 < 0 (after background subtraction) for oscillation 
studies. For each topology (see Fig. 1 on the left) the number of measured events, the MC prediction for 
no-oscillations and the ratio (Data/MC-no osc) are given §, i, §• 


I thank the members of the MACRO collaboration, in particular G. Giacomelli and M. Spurio. 
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